Control of coherence resonance in semiconductor superlattices 
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We study the effect of time-delayed feedback control and Gaussian white noise on the spatio- 
temporal charge dynamics in a semiconductor superlattice. The system is prepared in a regime 
where the deterministic dynamics is close to a global bifurcation, namely a saddle-node bifurcation 
on a limit cycle (SNIPER). In the absence of control, noise can induce electron charge front motion 
through the entire device, and coherence resonance is observed. We show that with appropriate 
selection of the time-delayed feedback parameters the effect of coherence resonance can either be 
enhanced or destroyed, and the coherence of stochastic domain motion at low noise intensity is dra- 
matically increased. Additionally, the purely delay-induced dynamics in the system is investigated, 
and a homoclinic bifurcation of a limit cycle is found. 



I. INTRODUCTION 



In contrast to the problem of controlling determinis- 
tic chaos, for which a number of methods have been 
proposed and successfully applied [l|, control of noise- 
induced and noise-mediated motion is a significantly less 
studied concept. Noise-induced dynamics refers to the 
case where there are no self-sustained oscillations in the 
absence of noise. The deterministic system rests in a 
stable steady state, e. g. a stable fixed point, and may be 
pushed out from it by random fluctuations. On the other 
hand, in the case of noise-mediated motion, the system 
already exhibits deterministic oscillations. Generally, the 
addition of noise may not only smear-out those oscilla- 
tions but may also induce qualitatively new structures 
and dynamics, like coherence resonance 0, S]- 

The ability to control the properties of noisy oscilla- 
tions (both noise-induced and noise-mediated) is very of- 
ten of practical relevance. This usually implies the en- 
hancement in the regularity of motion. However, in some 
cases, for instance in medical or biological applications, 
too much regularity is unwanted since it may lead to dam- 
aging consequences, e. g. epilepsy or Parkinson's disease 



Previous works have mainly concentrated on the con- 
trol of stochastic oscillations in low-dimensional sim- 
ple models [1 i, 0, i, i, E [nl, d, while control of 
noise-induced dynamics in spatially extended systems 
[3 [3 [3] seems still to be an open problem. This is the 
subject of the present paper. It is organized as follows. 
In Sect. II we will introduce the stochastic model of a 
semiconductor superlattice. In Sect. Ill time-delayed 
feedback control is implemented. In Sect. IV delay- 
induced spatio-temporal charge dynamics is studied in 
the absence of noise. Sect. V deals with the time-delayed 
feedback control of stochastic spatio-temporal dynamics. 



II. SEQUENTIAL TUNNELING MODEL 

We consider a stochastic model for a superlattice 
which consists of epitaxial layers of two alternating semi- 
conductor materials with different bandgaps, thus form- 
ing a periodic sequence of potential wells and barriers. 
The superlattice is a prominent example of a semicon- 
ductor nanostructure device which may serve as a prac- 
tically relevant nonlinear model system [l7j and may find 
applications as an ultra-high frequency electronic oscilla- 
tor [ll,[ll[2|. 

Our model is based on sequential tunneling of elec- 
trons plj . The resulting tunneling current density 
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^+l) from well m to well m + 1 de- 



pends only on the electric field between both wells 
and the electron densities n,„ and rim+i in the respective 



wells (in units of 



Ref. [21 1 



'), as given by eqs. (83), (86) in 



We simulate numerically a superlattice [24I of = 100 
GaAs wells of width w = 8nm, and Alo.3Gao.7As barriers 
of width b — 5nm, doping density Njj — lO^^cm"^ and 
at temperature T 20K 23] . 

The random fluctuations of the current densities are 
modelled by additive Gaussian white noise £,m(t) with 

iCmit)) - 0, {Umm'it')) = S{t - t')6mm', (1) 

which yields the following Langevin equations (to = 
l,...,iV): 



drin 



dt 



•.+D^rn{t) — Jm^m+l ^ D^„i+l{t), (2) 



where D is the noise intensity [l6|. Since the wells in 
our superlattice model are considered to be weakly cou- 
pled and the tunneling times are much smaller than the 
characteristic time scale of the global current 
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these noise sources can be treated as uncorrelated both 
in time and space. Charge conservation is automatically 



guaranteed by adding a noise term to each current 
density Jm-i^m- 

The electron densities and the electric fields are cou- 
pled by the discrete Poisson equation: 

ereo(^m - F„i-i) = e(n,„ - No) for m = 1, ... TV, (4) 

where and eq a-re the relative and absolute permittiv- 
ities, e < is the electron charge and Fq and Fjv are 
the fields at the emitter and collector barrier, respec- 
tively. In the deterministic case, i. e. a,t D — 0, rich dy- 
namical scenarios can be observed including formation 
of charge accumulation and depletion fronts associated 
with field domains bounded by these charge fronts, and 
with stationary, periodic or even chaotic current oscilla- 
tions [22, |2J, |25| . The two control parameters which are 
crucial for the dynamics, are the applied voltage U be- 
tween emitter and collector, which gives rise to a global 
constraint: 

N 
m=0 

where d is the superlattice period, and the contact con- 
ductivity cr For the current densities at the emit- 
ter and collector we use Ohmic boundary conditions, 
Jo^i = (tFo and Jn^n+i = aF^jf^- 

III. TIME-DELAYED FEEDBACK 

Time-delayed feedback was originally proposed by 
Pyragas in the context of chaos control [27|- The goal 
of this method is to stabilize unstable periodic orbits 
(UPOs) by adding, to a chaotic system, a control force 
in the form of the difference between a system variable 
at time t and its delayed version at time t — t. Here t is 
a delay time and K is the feedback strength. 

Apart from the control of UPOs the stabilization of un- 
stable fix ed p oints can also be achieved by time-delayed 
feedback [28l. [29| . This method proved to be very power- 
ful and has been successfully applied to various physical 
systems since then Other variants have been elab- 
orated, e. g. extended time-delay autosynchronization 
(ETDAS) and have been applied not only to deter- 
ministic systems (sil . [s^ Issl . [SJ] including fixed points 
but to stochastic systems as well [s^ . 

Apart from the deliberate design of feedback control 
loops to manipulate the intrinsic dynamics, delay arises 
also naturally in many systems. Typical examples of such 
systems are lasers where the delay enters through the 
coupling to external cavities (optical feedback) [sj [H, 
Isgj and neurons, where the signal propagation yields a 
delay time [40| . 

An easy way to implement control in the superlattice 
model is to choose the output signal to be the total cur- 
rent density, J[t) and simply add the control force to the 
external voltage U, i. e. 

U = Uo-K(j{t)-Jit-T)), (6) 
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FIG. 1: (color online) Schematic phase portrait below a 
SNIPER bifurcation. The full circle and the open circle mark 
the stable node and the saddle-point, respectively. 

where Uq is the time-independent external voltage bias 
and J is the total current density, filtered by a low-pass 
filter, 

J{t) f J(i')e~"**"*'^di', (7) 
Jo 

with a being the cutoff frequency ^33j . 

Since both voltage and total current density are exter- 
nally accessible global variables, such a control scheme 
is easy to implement experimentally. The low-pass filter 
in the current density was introduced initially for the ef- 
fective control of chaotic motion in the superlattice [s^ . 
Due to the well-to- well hopping of the depletion and ac- 
cumulation fronts, the current is modulated by fast small- 
amplitude oscillations. These high-frequency oscillations 
render the control loop unstable as J{t) is fed back to 
the system and therefore they need to be filtered out. 
Stochastic oscillations, like chaotic ones, exhibit the same 
high-frequency oscillations and thus we apply this filter 
here, too. 



IV. DELAY-INDUCED DYNAMICS 

In the following, the system is prepared in a stable 
fixed point corresponding to values of voltage U = 2.99U 
and conductivity a = 2.082J7~^m^"'^. It corresponds to 
a stationary accumulation front separating the high-field 
domain attached to the collector from the low-field do- 
main adjacant to the emitter. Besides the stable fixed 
point (node) there exists a saddle-point, whose unstable 
manifolds are connected to the stable node, forming a 
closed loop in phase space as depicted schematically in 
FigH] [16]. In this regime the system is very close to a 
saddle-node bifurcation on a limit cycle (saddle-node in- 
finite period bifurcation, SNIPER). In [l^ a bifurcation 
analysis in the (a, U) plane was performed showing how 
SNIPER bifurcations govern the transition from station- 
ary to moving field domains in the superlattice. In the 
vicinity of such a bifurcation the system is excitable [4l| 
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and therefore very sensitive to noise which is able to trig- 
ger front motion through the device. Moreover, the phe- 
nomenon of coherence resonance Q was also confirmed 
in our model. 

In this paper, we will study how time-delayed feedback 
acts on the system both in the presence and absence of 
noise, when prepared in the vicinity of a SNIPER bifur- 
cation. 

In the absence of both delay and noise the only stable 
attractor in phase space is the stable node. Regardless 
of the initial condition, all trajectories end there. In par- 
ticular, when selecting the initial condition on one of the 
two unstable manifolds of the saddle-point, the system 
performs a long excursion along the invariant circle be- 
fore ending in the stable node (Fig[T|). This deterministic 
path is affected by the delay for given combinations of the 
two control parameters, K and r. By keeping the delay 
T fixed and considering different control amplitudes K 
we track these orbits in phase space (Fig. [2]). The top 
panel of Fig.[2lla) (K=0) shows a trajectory which closely 
follows the unstable manifold of the saddle-point (cross) 
and ends in the stable node (full circle), cf. blow-up in 
the middle panel. Note that in the chosen 2-dimensional 
projection of the N-dimcnsional phase space the invari- 
ant circle (saddle-node loop) is distorted to a figure-eight 
shape. 

In Fig.[2)^b), during the first few nanoseconds the sys- 
tems acts as it would in the absence of delay, repelled by 
the saddle-point. Control is switched on at t = t = 2ns 
when the control force begins to act. The interval [0, 2ns] 
serves as initial condition of the delay equation. This 
becomes evident in a "twist" in the trajectory just be- 
fore the orbit reaches the stable node (middle panel of 
Fig. El^b)). For a moment it looks as if the system is at- 
tracted to the saddle-point instead of the node. This may 
be understood as follows: The control force shortly pulls 
the system off the phase space of the uncontrolled sys- 
tem pushing it towards the stable manifold of the saddle- 
point. At a critical value = 0.00QA375Vmm'^ /A, the 
system is indeed "swallowed" by the saddle and the tra- 
jectory closes in a homoclinic orbit. 

In the top panel of Figs. [2Kb) - ^d) the trajecto- 
ries for three values of K approaching this critical value 
are shown. Due to the high dimensionality of the sys- 
tem, which is blown up to infinity due to the delay, the 
above mechanism is not clearly demonstrated in a 2- 
dimensional projection in phase space. One must zoom 
in carefully in order to see the deviation from the de- 
terministic path due to delay (middle panels of Fig. 
This deviation is even better visible in the bottom pan- 
els of Fig. [2] where the final part of the electron density 
time series of ngs is plotted. In (a) the deterministic 
trajectory is plotted and the thick grey solid and dashed 
lines mark the position of the stable node and the saddle- 
point, respectively. It is clear that, the closer one is to 
the homoclinic bifurcation (d), the closer to the saddle- 
point does the system reach and the longer the trajectory 
remains there, before ending up in the stable node. 




FIG. 2: (color online) Top: Phase portrait in the (n65,n64) 
plane for fixed r = 2ns and (a) K = 0, (b) 0.006, (c) 0.0064 
and (d) 0.006437V mm^ /A. Middle: Enlargement of corre- 
sponding phase portraits shown in upper plots. Cross marks 
the position of the saddle-point, full circle marks the posi- 
tion of the stable node and IC is the initial condition, which 
is the same in all figures. Bottom: Corresponding electron 
density time series (only the part approaching asymptotically 
the stable node of nes). Thick grey line denotes the value 
of the stable node and dashed grey line denotes that of the 
saddle-point. U = 2.99V and a = 2.0821(r2m)"\ 



Beyond the critical value of Kc — O.OOGMVmm'^ /A 
the homoclinic orbit breaks and a limit cycle is born. 
In Figs. [3Ka) - [3Kc) the time series of the current den- 
sity for three different values of K above Kc are shown. 
The period clearly decreases for increasing K. Plotting 
it as a function of the control strength we obtain Fig. [3] 
(d) . The period T shows a characteristic scaling law 
|42j . T ~ ln[K — Kc], shown in the inset. This law gov- 
erns another type of global bifurcation, namely the ho- 
moclinic bifurcation. This delay-induced dynamics is in 
perfect agreement with a generic two- variable model for 
the SNIPER bifurcation, which also exhibits homoclinic 
bifurations of a limit cycle if delay is added [i^. Here, 
however, we have a much more complex spatio-temporal 
system. The corresponding space-time plots are shown in 
Fig. [31 It is clear that accumulation and depletion fronts 
corresponding to dipole domains are created at the emit- 
ter (bottom), and move through the device. Due to the 
global voltage constraint, they interact with the addi- 
tional accumulation front in the middle of the sample, 
thus forming a tripolc oscillation [22]. 

Following the homoclinic bifurcation in the K-t plane 
we numerically obtain the regime where control induces 
limit cycle oscillations. The result is shown in Fig. [5] 
Crossing the bifurcation line by increasing either K or r, 
oscillations are born, whose period decreases with K or 

T. 
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FIG. 3: (a) Time series of the delay-induced limit cycle for 
fixed T — 2ns showing period lengthening as the homo- 
clinic bifurcation is approached from above {K > Kc), for 
K = 0.0064379, (b) K = 0.007 and (c) K = O.OlWmm'^ /A. 
(d) Period T as a function of control amplitude K, and char- 
acteristic scaling law governing the homoclinic bifurcation 
0.00643751/mmV^, U = 2.99^^^, and 



shown in inset. 
a = 2.0821(r2m) 
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FIG. 4: (color online). Delay-induced front motion; Space- 
time plots of the electron density for K = 0.0064379, (b) 
K = 0.007 and (c) K = 0.019FmmV^. Light (red) and 
dark (blue) shading corresponds to electron accumulation and 
depletion fronts, respectively. The emitter is at the bottom, 
the collector at the top. Parameters as in Fig. |31 



CONTROL OF COHERENCE RESONANCE 



In this section, Gaussian white noise is added accord- 
ing to Eqs. H]), By considering two values of the 
control strength which lie inside or outside of the regime 
where delay induces a limit cycle (FigO), we will study 
the effect of the time delay on noise-mediated and noise- 
induced oscillations, respectively. 

The regularity or coherence of noisy oscillations can 
be quantified by various measures. Here we use (i) the 
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FIG. 5: (color online) Bifurcation diagram in the K-t plane 
below the SNIPER bifurcation in the superlattice. Yellow 
(light) area denotes the parameter regime for which delay- 
induced oscillation occur. In the green (dark) area there 
are no delay-induced bifurcations. The black line marks the 
delay-induced homoclinic bifurcation line. U = 2.99^ and 
a = 2.0821(nm)-^ 



correlation time 44] : 



ir.nr — 



1 



^-(0) Jo 



\ip{s)\ds, 



(8) 



where 7p{s) = ({J{t) - (J))(J(t + s) - (J))) is the au- 
tocorrelation function of the current density signal J(t), 
and (ii) the normalized fluctuation of pulse durations [3| : 



Rt — 



\J{T^)-{Ty 

(T) 



(9) 



typically used for excitable systems exhibiting oscilla- 
tions in the form of spike trains with two distinct time 
scales. These time scales are the activation time, which is 
the time needed to excite the system from the stable fixed 
point and the excursion time which is the time needed to 
return from the excited state to the fixed point. The sum 
of these two times equals the pulse duration or period of 
the oscillation T, which denotes the time between two 
spikes and is also known as interspike interval. 

By keeping the noise intensity fixed ai D = 
l.OAs^/^/m^ we first select a value K = 0.002Vmm'^/A 
outside the delay-induced limit cycle regime. This cor- 
responds to the lower dashed horizontal line in Fig. O 
In the right panel of Fig. Oa), the correlation time is 
plotted versus the time delay. It exhibits a slight mod- 
ulation with a period close to the period of the noise- 
induced oscillations {T) — 14.5ns and reaches min- 
imum values for r = n(T),n e N. Overall, however, it 
remains close to the control-free value, tcor = 19.76ns. 
At if = 0.02Vmm^ /A inside the delay-induced limit cy- 
cle regime (upper dashed line in Fig. , this modulation 
is much stronger and has a period close to the delay- 
induced period (T = 11ns, see Fig. [Sl[d)). In addition. 
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FIG. 6: (color online) Mean interspike interval (T) (left) 
and correlation time tcor (right) in dependence on the time 
delay r. (a) Control strength K — 0.002Vmm^ /A and 
noise intensity D = l.OAs^^'^ /m'^ , (b) K = 0.02Vmm'^/A 
and D = l.OAs^^^/m'^ and (c) K = 0.02Vmm'^/A and 
D = 2.5As^^'^/m^. Averages over 30 time series realizations 
of length T = 1600ns have been used for the calculation of 
averages over 1000 periods for (T). U — 2.99V and 
a = 2.0821(nm)-^ 
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one can better distinguish between non-optimal and op- 
timal values of r at which the correlation time attains 
maximum values. This is shown in the right panel of 
Fig. [6] (b). For a higher noise intensity (Fig. [6] (c), right 
panel) the effect is similar but weaker. 

Next we are interested in how the time scales are af- 
fected by the delay. We express the time scales through 
the mean interspike interval (T) and look at its depen- 
dence upon the time delay r for a fixed value of the 
noise intensity, D = l.OAs^^^ /m^ , and control strength 
K = Q.002Vmm^ /A, chosen outside of the delay-induced 
oscillations regime (lower dashed line in Fig. [5]). As 
shown in the left panel of Fig. [5] (a), (T) is slightly mod- 
ulated due to the delay with a period close to the noise- 
induced mean period ((T) « 14.5ns) [l6j . 

In the left panel of Fig. [5] (b) a value of K in- 
side the delay-induced oscillations regime is used, K = 
0.02Vmm^ /A (upper dashed line in Fig. For t — 0, 
the mean interspike interval is equal to the noise-induced 
period, (T) 14.5ns [l^. As the time delay increases, 
and the delay-induced bifurcation line is crossed, (T) 
sharply drops to the value of 11ns which corresponds 
to the period induced by the delay (see Fig. [Hd)). By 
further increase of r, (T) rises a little above 12ns. Then, 
for T = 11ns the mean interspike interval decreases again 
and the same scenario is repeated with a modulation pe- 
riod very close to the delay-induced period. There is 
some resemblance to the piecewise linear dependence of 
(T) upon T reported in other excitable systems: The 
FitzHugh-Nagumo model in d, 0, US] and the Oregona- 
tor model of the Belousov-Zhabotinsky reaction (under 
correlated noise and nonlinear delayed feedback) in [l^ 
which, like our system, is also spatially extended. The 
difference to our present analysis is that in those mod- 



FIG. 7: (color online) Correlation time (right) and normalized 
fluctuation of pulse durations (left) as a function of the noise 
intensity for (a) K = 0, (b) {K,r) = (0.02T/mm^/A, 11ns) 
and (c) {K,t) = (0.02ymmV^, 14.5ns). All three curves 
plotted together in (d). Averages over 30 time series realiza- 
tions of length T = 1600ns have been used for the calculation 
of tcor and averages over 1000 periods for Rt- U — 2.99V^ 
and a = 2.0821(r2m)"\ 



els the case of a delay-induced limit cycle was excluded. 
An explanation for the entrainment of the time scales by 
the delayed feedback in case of systems below a Hopf bi- 
furcation d, 0, H, [l3| was given on the basis of a linear 
stability analysis. It was shown that the basic period is 
proportional to the inverse of the imaginary part of the 
eigenvalue of the fixed point which itself depends linearly 
upon r, for large time delays. The effect of noise and de- 
lay in excitable systems was also studied analytically in 
pJOl . IT3 | based on waiting time distributions and renewal 
theory. 

Finally we look into the dependence of the correlation 
time tcor and the normalized fluctuation of the interspike 
intervals Rt on the noise intensity. We keep the control 
strength fixed to the value corresponding to Fig. [SJb) 
(right panel), from which we also select an optimal and 
a non-optimal value of the time delay and compare the 
results to the uncontrolled system. In the left and right 
panel of Fig. [71 Rt and tcor are plotted, respectively. 
The case K — is shown in Fig. [TJa) for direct com- 
parison. Coherence resonance shows up as a minimum 
of Rt and a maximum of tcor, respectively. For both 
non-optimal (Fig. [7]^b)) and optimal r (Fig. [Tfc)), there 
is an enhancement in the coherence at low noise inten- 
sity. Correlation times attain much higher values than 
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those of the uncontrolled system, especially at low noise 
level. Similarly, the interspike interval fluctuation is 
much smaller. In addition, for non-optimal delay time 
T « 11ns, the effect of coherence resonance is suppressed 
(Fig. [Ztb)). The correlation time still shows a small local 
maximum exactly where the uncontrolled system does, 
but for small noise intensities the correlation time dra- 
matically increases in a monotonic way to much larger 
values of tcor- On the other hand, for optimal r k, 14.5ns 
(Fig. Efc)), coherence resonance is maintained and both 
tcor and Rt show a maximum and minimum, respec- 
tively, but at a much lower noise intensity than in the 
free system. The comparison between all three cases is 
better visible in Fig.[71^d) where the three curves are plot- 
ted together. 

VI. CONCLUSIONS 

We have shown that by applying a time-delayed feed- 
back force to a semiconductor superlattice stationary 
field domains (bounded by charge accumulation and de- 



pletion fronts) can be transformed into travelling do- 
mains in a homoclinic bifurcation of a limit cycle if the 
system is prepared below a saddle-node infinite period 
SNIPER bifurcation. With the addition of Gaussian 
white noise, control results in a modulation of both the 
coherence and time scales of the system with the time 
delay. The periodicity of this modulation is determined 
by the competition between the different time scales im- 
posed by noise and control and their dependence on the 
noise intensity and time-delay, respectively. We distin- 
guish between optimal and non-optimal time delays at 
which the coherenc resonance effect is enhanced or de- 
stroyed, respectively. In both cases the correlation times 
of stochastic domain motion are dramatically increased 
at low noise intensities. 
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